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Several features of an f(R) theory in which there is a maximum value for the curvature are 
analyzed. The theory admits all the vacuum solutions of General Relativity, and also the radiation 
evolution for the scale factor of the standard cosmological model. Working in the Jordan frame, a 
complete analysis of the phase space is performed, and shown to be in agreement with examples 
obtained by numerical integration. In particular, nonsingular cosmological solutions were obtained, 
which after the bounce enter a phase of de Sitter expansion and subsequently relax to a GR-like 
radiation-dominated evolution. 
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I. INTRODUCTION 

Although General Relativity (GR) agrees with observation with good precision in several situations pQ, theories 
described by a Lagrangian that is a function of the Ricci and other scalars constructed with the Riemann tensor have 
been intensively studied in recent times. These theories, which are carefully engineered in such a way that they differ 
from GR either in the limit of low or high curvature (or in both limits) have a long tradition, starting with a paper 
by Weyl in 1918 0. The current interest in these theories is twofold. First -in the case of low curvature- the aim is 
to model the accelerated expansion of the universe that follows from several observations [3J (when interpreted in the 
standard cosmological model [1]) without using dark energy [5]. Second, since we have yet no direct observational 
evidence of the behaviour of the gravitational field for high values of the curvature, these theories are used to widen 
our understanding of possible new features of gravity in this regime. Black holes [B] and neutron stars [7] are natural 
candidates to look for deviations from GR in the strong regime, as well as the universe in its phase of very high 
density. We shall be interested here precisely in the latter circumstance, by restricting to theories with Lagrangian 
that is a function of the Ricci scalar only [3]. This option is due to simplicity, and to the fact that f(R) theories are 
favoured over more general choices because they do not suffer from the Ostrogradski instability [5]. 

The cosmology of f(R) theories in the large-curvature limit has been explored for instance in |10j . where inflation 
was analyzed in the theory f(R) = R + R 2 /M 2 . Still in this regime, f(R) theories have been used to remedy one of 
the main problems of the standard cosmological problem, namely the initial singularity |11) . In this regard, general 
conditions for the absence of an initial singularity in f(R) theories were discussed in [12], while solutions displaying 
a bounce for specific forms of f(R) were exhibited in [15] . 
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In this article, we shall begin the analysis of the large-curvature characteristics of the theory given by 



/ R 2 

f(R) = R-0jl--p+P, (1) 

where /3 is a constant that explicitly sets un upper limit for the curvature, in a way reminiscent of that of Born-Infeld's 
electromagnetism [14]. Notice however that our proposal is different from the so-called Born-Infeld gravity, which 
has been discussed for instance in |15) . The theory examined here is also different from the one studied in |16j . with 
Lagrangian given by 

C =^ ^_ (2) 

16vrG y/\ + WW ' K} 

(where I is a dimensional constant) since in this theory the curvature scalar is not bounded and the effective gravi- 
tational constant G c ff (which is proportional to ^) diverges for large R, while in the theory defined by Eqn.|l| the 
curvature is limited by the parameter /3, and G c g is finite for all R . 

Our main goal will be to characterize the cosmological solutions of the theory given by Eqn.Q in the strong field 
regime (in the presence of radiation), giving particular attention to solutions that do not display an initial singularity. 
The general features of this theory, as well as its restriction to the standard cosmological model will be discussed 
in Section [llj A phase space analysis of the resulting equations will be presented in Section III [17] , while relevant 



examples of this analysis will be illustrated by numerical integration in Section IV We close with a discussion of the 
results in Section IVl 



II. GENERAL FEATURES 



The starting point is the gravitational action given by 

1 

2k 

where k = 8ttG, G being the gravitational constant, we have set c = 1, and f(R) is a general function of the Ricci 
scalar. As mentioned above, we shall work with the function 



S g = — I d*xV=gf(R), (3) 



/ R 2 

f{R) = R-pyl\--p+p. (4) 

The form above is similar to the lagrangian proposed by Born and Infeld to solve the problem of singularities in the 
electromagnetic field, and so it carries their name. /3 is a free parameter of the theory, hopefully to be identified with 
a fundamental constant of nature, like the inverse of the Planck length squared. This function bounds the value of the 
Ricci scalar, both from below and from above, with \R\ < /3, therefore singularities characterized by the divergence 
of R should be supressed in the solutions of this theory. The third term is inserted to recover Minskowski spacetime 
in the absence of curvature. The equations of motion that follow from the variation of the action ^ with respect to 
the metric are, 

^r*W ~ \f(R)%» - [V„V„ - g, v U] ffl = kT„ v , (5) 
where T^„ is the energy-momentum of the matter fields, defined by 

rp 2 SS M 

and the covariant derivative is defined using the usual Levi-Civita connection. Taking the trace, we obtain 

^-2/( fl ) + 3D«S=«r, (0) 
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We see from Eqns.([5]) and ^ in the case of the f(R) given by Eqn.Q that all solutions of GR with a traceless 
energy-momentum tensor as a source are solutions of the f(R) theory studied here, although they may not be unique. 
The same can be said about the vacuum solutions of GR. 

Regarding the two basic conditions that a given /(i?) should satisfy in order to be a viable theory, namely df jdR > 
and d 2 f/dR 2 > jSJ, a simple calculation shows that the second one is satisfied for any value of \R\ < f3, while the 
first one is violated for R < —0.71/3. However, considering that /3 is of the order of the inverse of Planck length 
squared, this violation takes place when the scalar curvature scale of spacetime, defined as l c = |i?| -1 / 2 , approachs 
the Planck length, l c > ^pianck, where quantum gravitational effects migth be important and the theory should be 
modified. We shall assume that the theory has plenty of solutions that do not reach this regime, and show below that 
this assumption is well-founded in a cosmological context. 

Since R contains second derivatives of the fundamental coordinate g^u , the metric tensor, the field equations arising 
from it are of fourth order in the derivatives of g M „. Both the analytical and numerical handling of the resulting 
equations is easier if we introduce the auxiliary field s = R, in such a way that the action is now written as 



where is a Lagrange multiplier. Variation of the action ^ with respect to s yields <f> = f'(s), where the prime 
denotes derivative with respect to the argument. The equation Q can be rewritten as 

Sg = hJ foV^alW-vtt)] > v(<t>) = <M0) - M4>)], (8) 

so the f(R) theory is identical to a Brans-Dicke theory with parameter cj = 0. The relation between 4> and s is as 
follows: 



The variation of ^ with respect to cf> results in 

R = V'{4>), (10) 
It is possible to see that, in order to be consistent with Eq. (fTol) one must set, 



a (4,) = p—J t 1 (11) 

Vi + (0-i) 2 

Hence, the potential must have the form 

v(<t>) = /3 (yr+ (0-1) 2 - 1) (12) 

In the situation of zero curvature, <j) = 1 and V(<j>) — 0, and Minkowski spacetime is regained. Finally, the variation 
of the action ^ with respect to the metric tensor. We obtain the following field equations 

4>G^ + -Vi^g^u - V^V,,!?!) + g^O<j) = kT^ v , (13) 



and the associated trace is 



□0-I^'(^) + ?^) = |T, (14) 



where we used Eqn.( 10 ). Matter has been included in the right hand side of these equations, as the energy-momentum 
tensor T^ v , with trace T. 

Armed with these equations, we turn next to their action on the geometry of the standard cosmological model. The 



background evolution is obtained by inserting the Friedman-Lemaitre-Robertson- Walker (FLRW) metric in Eqns.(13) 
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and (14 1, and considering a perfect fluid form for the energy-momentum tensor. Then, we get for the 00 and 11 

V p 



components and for the trace, respectively, 



H z + H<P = 



6 



(2H + 3H 2 



V 

~-2H<t>- 



P 



(15) 
(16) 
(17) 



Since we are interested in the early universe, which was radiation dominated, the right-hand-side of the trace equation 
is zero, and we can write it as the Klein-Gordon equation for the field <p 



4> + 3H<p + U'(<f)) =0, 



with the effective potential 

U(<t>) 



</Vl + (</>- l) 2 - 3-v/l + ((/)- l) 2 + 3 arcsinh((/> - 1) - 40 



(18) 



(19) 



shown in Fig{T] This plot justifies a posteriori the use of the so-called Jordan frame, since the potential is not 
multivalued. 




FIG. 1: The effective potential U(<f>) displays a local minimum at <f> = 1, corresponding to GR, and a local maximum 

at 6 = 3.315. 



The system of coupled equations (15l-(17l is very complicated to be solved analytically. The plot shown in Fig. [T] 
will be important to qualitatively understand the different behaviors that will be obtained when we numerically solve 
our system of equations. We will see that the local maximum and minimum in the potential correspond to critical 
points of the dynamical system associated the evolution equations. 



III. DYNAMICAL SYSTEM ANALYSIS 



The dynamical system is built from Eqns.( 15 )-( 17 1 and (|T8|) and is given by 

-2H 2 



H 

y = - 
<j> = y 



V 



6 ' 

3Hy - h$V' + -V, 



(20) 

(21) 
(22) 
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Here the curvature scalar R has been expressed in terms of the scale factor and its derivatives, taking the Hubble 
parameter as the fundamental function, and we introduced the variable y so that the Klein-Gordon equation can be 
reduced to two first order equations. 



A. Critical points 



The critical points, which correspond to those points where H = y 



H = ±y/V'/12, 4>V' = 2V, 
From the relation <j)V' = 2V, we find the following equation values for 



= 0, are given by 

y = o. 



u 3 - 2u 2 



= 1, 

it -4 = 0, 



(u = 0-l). 



(23) 

(24) 
(25) 



So we have in principle four values of 0, leading to 8 critical points if we consider H = ± \JV /12. But the third order 
polynomial equation has only one real root, and so we have in principle four critical points. Also, 0=1 corresponds 
to H = V = V' = 0. Hence we end up with three critical points given by: 



Pi 
Pi 
P 3 



= 1, H = 0, y = 0; 



= 3.315, H = + y / V , /12, y = Q; 
= 3.315, H = -tJV/12, y = 0. 



(26) 
(27) 
(28) 



The first point, Pi, corresponds to the Minkowski spacetime, while Pi to a expanding de Sitter and P3 to a contracting 
de Sitter, both with a constant scalar curvature R = 0.92/3. Notice that the values obtained for </> correspond to the 
local minimum and maximum in Fig. [T] To find the nature of these critical points we must linearize the system 



(|20|)-(|22]), by replacing 

H = H a +6H, = 0o + ^, V = Sy, (29) 
where Hq and 0o are the values of H and at the critical points (notice that yo — for all the critical points), given 



by Eqns.(26l-(28). We then find the following linearized system 

1 



SH 

8y 



-4H SH 



G 



-V"t 



1, 



-3HSy -3ySH - -V'i 
3 



-0F"50+ -V'i 



(30) 
(31) 

50 = Sy, (32) 
where we have omited the subindex on the background quantities. We shall investigate next each point separately. 
• Critical point Pi 

Corresponding to Minskowski spacetime, this point has Hq = and 0o = 1. The equations around this point 
are 

P CI 

-3*0, 
= Sy, 



SH 
Sy 



(33) 

(34) 
(35) 



which in matrix form read 









( SH 







Sy 


1 




\ 50 



(36) 
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The diagonalisation of this system involves the calculation of the eigenvalues A, given by 

/-AO § \ 
-A -§ 
1 -Ay 



det 



V 



= —A 3 — ^A = 



(37) 



The roots are 



A = ±i\ 



IP 
3' 



A = 0. 



(38) 



Since one eigenvalue is zero and the other two are pure imaginary with opposite signs, this critical point is a 
center-saddle point. 

Critical points P 2 and P 3 

These are de Sitter spacetimes. Following the same approach used for Pi, the linearized equations around the 
points P 2 and P3 are 



5H = -msH + -v"s<t> 

6 

Sy = -3H5y+^V'5cj)-^<j>V"54> 



Sip = Sy 

To obtain the eigenvalues, we must solve the equation 

1 



V + T} <t>V" A - -H(V' - 4>V") = 



X s + 7HX - 

For the critical point P 2 , we have 4> = 3.3146 and H = y/V'/12, therefore 

Ai = -1.10634, A 2 = -1.05456, A 3 = 0.224804. 



(39) 

(40) 
(41) 

(42) 
(43) 



Since we have two eigenvalues with negative sign and one with positive sign, this indicates that it is a saddle 
point in two different planes, while it is a attractor in another. 



For P 3 we put = 3.3146 and H = -y/V'/12, yielding 

Ai = -0.224804, A 2 = 1.05456, A 3 = 1.10634. 



(44) 



Once again, we have a saddle point in two different planes, but now it is an repeller in the third plane. We 
cannot tell which planes are those for now, but we will see that the repulsive/attractive behavior of these de 
Sitter points will be present in the numerical evolutions. 



B. Critical points at infinity 

The complete phase diagram can be determined from the knowledge of the critical points at infinity. Let us write 
a three-dimensional system as follows: 

x = X(x,y,z), (45) 

y = Y(x,y,z), (46) 

i = Z(x,y,z), (47) 



and define the new variables 



x= k> y= k> z = v (48) 
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with the constraint 

u 2 + v 2 + w 2 + k 2 = 1 (49) 

The new variables define a three-sphere, with k — denoting the "equator" of the three-sphere, in fact a two-sphere. 
The two-sphere defined in this way corresponds to the infinity. Hence, we have now 

y - %u = X{u/k,v/k,w/k) (50) 
k k z 

^ - ^v = Y{u/k,v/k,w/k) (51) 

y-^w = Z(u/k,v/k,w/k) (52) 
k k z 

with the condition 

udu + vdv + w dw + k dk = (53) 

Combining these equations we obtain the following system: 

u = {l-u 2 )kX -kuwZ -kuvY (54) 
v = (l-v 2 )kY -kvwZ -kuvX (55) 
w = (1 — w 2 )kZ — kuwX — kvwY (56) 

In our specific case, putting H — > x, <fi — > z and keeping y as it is, we have the expressions 

X- 1 / 2u 2 + P k 2 {w-k) \ 1 

*-fc*\ _2 " + 6 ^ + ( W -fc)^ '^ (5?) 

F = 1 ,y m - fc L + ^ [ V fc 2 + (™ - fc) 2 - ll 1 = ^ (58) 



k 2 1 3 ^fc 2 + - fc) 2 3 



fc 2 " 



z =tf*>=± i Z (59) 

Using the expressions for our specific system, redefining the evolutive parameter, such that dr/k — > dt, we end up 
with the following dynamical system at infinity (fc = 0): 

it = -2u 2 (l-u 2 )+3u 2 v 2 (60) 
v = -3uv{l-v 2 ) + 2u 3 v (61) 
w = 2u 3 w + 3uv 2 w (62) 

There are the following critical points: 

• P 4 : u = v = w = 0. Incompatible with the condition (at infinity) u 2 + v 2 + w 2 = 1. 

• P 5 : u = ±1, u = 0, w = 0. This critical point is located on the top and on the bottom of the u axis, that is, the 
axis of H at infinity, with the scalar field and its deriviative equal to zero. 

• Pq: u = 0, for all v, w satisfying the condition v 2 + w 2 = 1, or <p 2 + cf) 2 — > oo. 



The critical point P5 represents the beginning of the universe from a singularity (the big bang), sign plus, or the 
end of the universe in a singularity (big crunch), sign minus, the time reversal of the first case. 

The critical points Pq lying on the sphere v 2 + w 2 = 1 , or cb 2 + (fx 2 — > 00, correspond to the physical situation 



where <fi is deep down in the potential U((f>) given in Eq. (19). Looking at Eq. (20), this corresponds to contracting 



or expanding de Sitter spacetimes, where the Ricci scalar is saturated to its maximum constant value, R = f3 (that is 
the reason for having u = kH = k\/YlR = when k = 0). 

Linearizing the system around the critical point P5, we obtain, for that 

Sii = ±4Su. (63) 

This implies that the critical point u = 1 is a repeller (the curve emerges from this point) and u = — 1 is an attractor: 
the first critical point is the big bang while the second is the big crunch. 

Linearization of the equations around the critical point Pq, leads to the following system: 

5u = 0, (64) 
Sv = -3(l-v 2 )5u, (65) 
5w — 3v 2 w8u. (66) 

The first of these equations imply that 5u remains constant. Using the fact that u 2 + w 2 = 1, we obtain the system 

Sv — —3vw 2 Su (67) 



Sw = 3v w 6u (68) 



This implies the following configurations: 



• 8u > 
























- v > 0, 


w 


> 





implies 


that 


the trajectories 


are 


repulsive in the direction w and attractive 


in 


the direction 


v; 


- v > 0, 


IV 


< 





implies 


that 


the trajectories 


are 


attractive in both the directions v and u>; 








- v <0, 


w 


> 





implies 


that 


the trajectories 


are 


repulsive in both the directions v and w; 








- v < 0, 


w 


< 





implies 


that 


the trajectories 


are 


attractive in the direction w and repulsive 


in 


the direction 


v. 


• Su < 0: 
























- v > 0, 


IV 


> 





implies 


that 


the trajectories 


are 


attractive in the direction w and repulsive 


in 


the direction 


v; 


- v >0, 


w 


< 





implies 


that 


the trajectories 


are 


repulsive in both the directions v and w; 








- v < 0, 


w 


> 





implies 


that 


the trajectories 


are 


attractive in both the directions v and w; 








- v < 0, 


IV 


< 





implies 


that 


the trajectories 


are 


repulsive in the direction w and attractive 


in 


the direction 


v. 



After this analysis, one can envisage many scenarios. The less interesting one is the usual radiation-dominated, 
spatially-flat FLRW solution, which is allowed in this framework, coming from the critical point P 5 , u — 1 (the big 
bang singularity) and going to the critical point Pi, H = (Minkowski), always with cf> = 1. However, one may have 
many other singular solutions, e.g., one coming from the singular point, P5, u = 1, but away from the line (f> = 1 at 
infinity, which may cross the plane defined by H = 0, and go to the critical point at (f> — > —00 by below. 

The region that may exhibit nonsingular solutions is the one with positive (f>. Some of them can begin at the critical 
point Pi , contracting from Minkowski spacetime to the critical point P3 (contracting de Sitter spacetime) , then going 
to the critical point Pi (expanding de Sitter spacetime), with a bounce in between, where the plane H = is crossed. 
The trajectory then returns to the point Pi, that is, expanding to Minkowski spacetime. The scalar field starts to 
oscillate very close to 4> = 1 at t — > —00 (the general relativistic limit), goes to the local maximum of the potential (due 
to anti-friction because of contraction, see Fig{T]), and then returns to the local minimum at 4> — 1 at t — > —00. This 
is a nonsingular solution with a bounce and an almost de Sitter expansion happening before a decelerated expansion 
phase, close to the usual radiation-dominated, spatially-flat FLRW general relativistic solution. 

Other solutions are those contracting from the singular point Pq, a de Sitter spacetime with R as j3 and 4> — > 00 
at t — > —00, in which the scalar field climbs the potential from below, due to anti- friction. The solution crosses the 
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plane H = 0, thus going through a bounce, and reaches the critical P27 near the local maximum of the potential (see 
FigjTJ) , and then moves to the critical point P\, with « 1 at t — > 00. Again, one has a non singular solution with 
a bounce and an almost de Sitter expansion happening before a deccelerated expansion phase close to the radiation- 
dominated, spatially-flat FLRW general relativistic solution. The difference is that the former solution corresponds to 
a nonsingular bouncing solution between two Minkowski spacetimes, while the second corresponds to a nonsingular 
bouncing solution between de Sitter contraction and Minkowski spacetime. Such trajectories are presented in the 
next section. 



IV. NUMERICAL RESULTS 



The dynamical system ( 20 )-( 22 1 can be numerically solved by simple inputs in programs like Maple or Mathematica. 
However, this system is highly nonlinear and thus very sensitive to initial conditions. We will begin our evolution 
around different critical points and see how the trajectories behave for small fluctuations (5H, 8<j> and Sy). For 
numerical considerations, we adopt (3 = 0.5. 

From Figfl] we see that the critical points of the finite region are the local minimum (Minkowski configuration) 
and maximum (de Sitter configuration). We will now list interesting examples obtained from particular sets of H, cf>, 
y and their fluctuations. 

• Around critical point Pi 

Starting from a Minkowski spacetime, the value of (f> begins to oscillate with increasing amplitude, while H 
oscillates around a decreasing average value. Eventually, cj> overcomes the local maximum at its right and rolls 
down the potential, towards (f> — > +00. This leads to a big crunch singularity, with H — > —00. The evolution is 
shown below. 




FIG. 2: Evolution for H and </> around Minkoski spacetime. The detailed frame shows the oscillation of the Hubble 

parameter between t = —80 to t = —20. 



• Around critical point P2 

Starting our evolution around the local maximum is very tricky, since fine tunning of 5<j) and Sy is needed 
in so that the field does not roll down the potential to 4> ~ * +00. An interesting solution is shown below 
(Figj3|, representing an universe emerging from a past singularity, as in the big bang model, and reaching the 
Minkowski configuration after an intermediate de Sitter phase, which is a saddle point. This is precisely the 
evolution presented in the cosmological Standard Model. 

• Around critical point P3 

In this case, we have the same situation as in P%, with the substitution H — > —H and t — > —t. The universe 
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FIG. 3: Evolution of a universe that emerges from a past singularity, enters a de Sitter phase and reaches 

Minkowski spacetime. 



begins in a Minkowski state, but the Hubble parameter begins to oscillate. Then, it goes through a contracting 
de Sitter phase which culminates in a big crunch singularity, with H — > — oo 




FIG. 4: In this case the Hubble parameter oscillates around the Minkowski state, and goes through a contracting de 

Sitter phase which culminates in a big crunch singularity. 



• Bouncing solutions 

— Minkowski-to-Minkowski evolution: taking the starting point of our evolution as the intermediate point 
between both de Sitter configurations, with <f> ss 3.315 and H we obtain trajectories characterizing 
nonsingular universes. Two examples are shown below. The first one begins in a Minkowski state which 
oscillates with growing amplitude, until it reaches a contracting de Sitter phase. Suddenly, the universe 
bounces to a de Sitter expansion, which lasts for a while. Then, it goes through oscillations once again 
to finally reach a Minkowski configuration (see Fig[5|. This is the only solution entirely contained in the 
finite region, with tf> always to the left of the local maximum of the effective potential shown in Fig[T] 

— de Sitter-to-Minkowski evolution: the second situation (Figj6]) also presents a bounce, but it starts from a 
contracting de Sitter state at t — > — oo. This reflects the result we obtained from the analysis of the critical 
point P3, which behaves as a repeller in this particular plane. 
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FIG. 5: The model in the figure begins with growing oscillations around the Minkowski state, followed by a 
contracting de Sitter phase, a bounce to a de Sitter expansion, to finally reach again the Minkowski geometry. 




FIG. 6: In this case the evolution, having started from de Sitter spacetime, displays a bounce, and reaches 

Minkowski in the end. 



— de Sitter-to-de Sitter evolution: once again, we get the repulsive and attractive behaviors of the de Sitter 
critical points in particular trajectories. The universe comes from an infinite de Sitter contraction in the 
past and goes through a bounce, acquiring afterwards a forever expanding de Sitter configuration (Fig[7]) . 
The field <fi in this evolution is restricted to the right of the local maximum of the effective potential, i.e., 
4> > 3.315 for all times. 

A bouncing solution can also be obtained analytically by assuming that the Ricci scalar approaches its maximum 
possible value, given by jH, at the bounce (t = 0), so we may write 

R = P-e, < e « 1. (69) 

Since e is very small, the component 00 of Eqn. (pi) in a FLRW universe is now given by 
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H(t) 




80-60-40-20 20 40 60 




FIG. 7: The figure shows a de Sitter to de Sitter evolution through a bounce. 



where we kept only the dominant terms, proportional to 1/Ve. Multiplying both sides by y2/3e, and assuming 
that the energy density at the bounce is high but finite, we get 



a 
a 



at 
2^~e 



oc 



(71) 



Expanding the scale factor in Taylor series around t — 0, and imposing the conditions a(t — 0) > and 
a(t = 0) = valid at the bounce, the form of the scale factor must be a(t) — a + a 2 t 2 + a 3 f 3 . Hence, the Ricci 
scalar is 



R = 6 



12<22 + SGa^t, 



(72) 



where we neglected orders higher than the linear, and we can identify j3 = 12a 2 and e = SGa^t, with a 3 < 0. 



However, inserting the scale factor in Eqn.(71) we obtain 

a he 



2ae 



= a 2 + 6a 3 t, 



(73) 



which is always different from zero and therefore our choice of a(t) cannot be a solution of the field equations. 
Although the presence of the term a^t 3 is inconsistent with equation (71), the same does not happen with a 



solution of the form a(t) = ao + a 2 t 2 + 04t 4 . The Ricci scalar in this case is 

R ~ 6 [2a 2 + (2«2 + 12a 4 )£ 2 ] , 

where we ignored terms of order higher than t 2 . It follows that /3 = 12a 2 and e 
we must impose 04 < — a 2 / 6. Inserting these quantities in Eqn.(71), we get 



a 
a 



ae 
2^ 



2a 2 



4a 2 t 2 
2t 2 



= 0, 



(74) 

-12(a| +6a 4 )t 2 . Since e > 0, 



(75) 



and so the scale factor expansion is indeed a solution. Extending the Taylor series to order t 5 , in the form 
a(t) = a + a 2 t 2 + a^ 4 + a 5 t 5 , it is straightforward to show that e = — 12(a 2 + Qa^t 2 + 120a 5 i 3 , yielding 



ae 
2ae 



cx t oc \/e, 



(76) 



and therefore ( 71 ) is satisfied. We can also see from the Taylor expansion that both a and a are finite at t = 0, 



so the Ricci scalar is well behaved in this scenario. We conclude that bouncing solutions can be obtained from 
the Lagrangian Q if the coefficient 03 in the Taylor expansion of the scale factor around t — is null. 
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V. FINAL REMARKS 

We have discussed here several features of an f(R) theory in which there is a maximum value for the curvature (in 
a way reminiscent of Born-Infeld electrodynamics) . The theory, defined by the Lagrangian given in Eqn. Q , admits 
all the vacuum solutions of GR, and also the radiation evolution for the scale factor of the standard cosmological 
model. Working in the Jordan representation, a complete analysis of the phase space was presented, the main results 
of which were illustrated with examples obtained by numerical integration. It was shown that there are cosmological 
solutions that attain a minimum value of the scale factor and then expand again, instead of displaying the so-called 
initial singularity. Perhaps the most interesting ones are those shown in Figs. ^ and ([6| in which after the bounce 
there is a phase of de Sitter expansion and a subsequent relaxation to GR, with the field <f> oscillating around = 1. 
In particular, the plots in Fig. ^ suggest that after the de Sitter expansion, the universe goes through a radiation- 
dominated evolution, at least for some time, before oscillating around the Minkowski stage. 1 

Although we worked in the Jordan frame, it is worth pointing out that we expect our main results to be valid also 
in the Einstein frame, whose metric is related to that of the Jordan frame via the transformation = ipg^. The 
relation between the Hubble parameter in the two frames is given by 

1 

Since H , <f> and its derivative (displayed in the plots) are well-behaved functions, it follows that the nonsingular models 
we obtained in the Jordan frame (all of which satisfy the condition <j) > for all values of t), will also be nonsingular in 
the Einstein frame. Moreover, since the conformal transformation is always regular, the instabilities due to transition 
from attractive and repulsive gravity |21j are absent in these bouncing models. 

Future work includes the study of the dynamical system with the addition of non-relativistic matter and dark 
energy [T5], and the investigation of the features of compact objects in this theory, such as neutron stars (along the 
lines presented in pI5]) and black holes. 
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